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Evolution of Weyl functions and 
initial-boundary value problems 

A.L. Sakhnovich 

Abstract 

This review is dedicated to some recent results on Weyl theory, 
inverse problems, evolution of the Weyl functions and applications 
to integrable wave equations in a semistrip and quarter-plane. For 
overdetermined initial-boundary value problems, we consider some ap¬ 
proaches, which help to reduce the number of the initial-boundary con¬ 
ditions. The interconnections between dynamical and spectral Dirac 
systems, between response and Weyl functions are studied as well. 
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1 Introduction 

Initial-boundary value problems for linear and nonlinear wave equations are 
interesting, difficult and have many applications since initial-bonndary con¬ 
ditions should be taken into acconnt in the stndy of various applied problems 
connected with wave processes. Weyl theory could be useful in this stndy 
(and, vice versa, initial-boundary valne problems are of interest in Weyl 
theory). A nnmber of related resnlts and references was presented in our 
book [58, Ch. 6]. Here we mainly review the results that appeared after the 
publication of [58] althongh some results from [58] and some earlier resnlts 
that are not contained in [58] are also discussed for the sake of completeness. 

In order to make the review reader-friendly and self-snfficient, we provide 
short schemes of the proofs for some of the main statements and start the 
review with the section ’’Preliminaries: Weyl fnnctions, inverse problems and 
zero curvatnre representation”. However, even this section contains some 
qnite recent important results. 

Zero curvature representation [3,26,44,63] is in a certain sense an ana¬ 
logue of the famous Lax pairs. It is a compatibility condition for linear 
systems, which are so called auxiliary systems for nonlinear integrable wave 
eqnations. This compatibility condition is discussed in Snbsection 2.3 of the 
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section ’’Preliminaries Here already, an initial-boundary value problem, 
which is basic for the study of the evolution of Weyl (Weyl-Titchmarsh) 
functions, appears. 

Section 3 is dedicated to generalized Weyl functions and corresponding 
inverse problems. In Section 4 and in Subsection 5.2 we consider some appli¬ 
cations of evolution formulas for Weyl (and generalized Weyl) functions to 
the problems of uniqueness and existence and to the unbounded solutions of 
wave equations. 

We note that Cauchy problems (in particular, Cauchy problems for in- 
tegrable nonlinear wave equations) are studied much more thoroughly than 
initial-boundary value problems. In fact, too many initial-boundary condi¬ 
tions are usually required for the study of initial-boundary value problems 
(see, e.g., [5,19,27,58] and references therein), and the problem becomes 
overdetermined. Thus, a crucial step here is the reduction of the initial¬ 
boundary conditions. Some important approaches to this reduction are dis¬ 
cussed in Section 5. 

Finally, Section 6 is dedicated to the initial-boundary value problem for 
dynamical Dirac system and its connections with Weyl theory. 

Appendix A contains some dehnitions and references from the theory of 
quasi-analytic functions. 

Notations. As usual, M stands for the real axis, M+ = (0, cxd), C stands 
for the complex plain, and C+ for the open upper semi-plane. We use no¬ 
tations Cm = {z : A( 2 :) > M} and C]b = {z : A( 2 :) < —M}. The equality 
D = diagjdi,...} means that D is a diagonal matrix with the entries di,... 
on the main diagonal and denotes the m x m identity matrix. ’’Locally” 
(e.g., locally summable) with respect to the semiaxis [0, oo) means: on all the 
hnite intervals [0, /] (e.g., summable on all [0, /]). The class of linear bounded 
operators acting from the normed space "Hi into the normed space 772 is de¬ 
noted by and we write simply BiTi) when "Hi = 772 =: 77. We 

say that a matrix function is boundedly differentiable when its derivative is 
bounded in the matrix norm, and an operator is boundedly invertible when 
the inverse operator exists and is bounded. The notation l.i.m. stands for 
the entrywise limit of a matrix function in the norm of L^(0,7), 0 < 7 < cxd. 
The notation DJI stands for the operator mapping a Weyl function (or gener- 
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alized Weyl function, depending on the context) (f(z) on the corresponding 
potential (e.g., = v for Dirac systems (2.1) and (3.1) where V has the 

form (2.2)). That is, is the solution of the inverse problem to recover 

the potential from the Weyl function. 

2 Preliminaries: Weyl functions, inverse prob¬ 
lems and zero curvature representation 

2.1 Main scheme 

One of the most important auxiliary linear systems in soliton theory is self- 

a; > 0, (2.1) 


mi + 1712 =■ (2.2) 

Irrii is the rrii x m* identity matrix and v{x) is an mi x m 2 matrix function. 

For instance, the well-known matrix defocusing nonlinear Schrodinger 
(defocusing NLS or dNLS) equation 

d 

2vt = i{v^x - 2vv*v) := —nj (2.3) 

admits [64,65] (i.e., is equivalent to) zero curvature representation 

Gt-Fx + [G, F] =0 ([G, F] := GF - FG), (2.4) 

G = i{zj+jV), F = -i{z^j + zjV-(iVx-jV^)/2), (2.5) 

where V has the form (2.2) and v = v{x, t). For this reason, systems Ux = Gy 
and yt = Fy are called auxiliary systems for dNLS, and it is easy to see that 
Vx = Gy has the form (2.1) for each hxed t. We note that V is called the 
potential of system (2.1). (Sometimes, for convenience, v is also called the 
potential.) 


adjoint Dirac (also called ZS or AKNS) system 

z) = i{zj + jVix))y{x, z), 

where 



0 

, V = 

■ 0 

V 

0 

I m 2 


V* 

0 _ 
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Weyl solutions are squarely summable on (0, cxd) solutions of Dirac system 
(2.1). Clearly, these solutions may be presented as linear combinations of the 
columns of the normalized by the condition 


n(0, z) = /, 


( 2 , 6 ) 


fundamental solution u of (2.1). Thus, Weyl-Titchmarsh (or simply Weyl) 
function of (2.1) is introduced in terms of u. 


Definition 2.1 Let Dirac system (2.1) on [0, cxd) he given and assume thatv 
is locally summable. Then Weyl function is an m 2 x mi holomorphic matrix 
function, which satisfies the inequality 


Imi ^iz)*]u{x,zyu{x,z) 


I \ dx < 00 , z E C+. 

Vpw) 


(2.7) 


The following proposition is proved in [29] (and in [58, Section 2.2]). 


Proposition 2.2 The Weyl function always exists and it is unique. More¬ 
over, the Weyl function is contractive, that is, ip{zY(p{z) < Imi- 


Inverse problems for the classical selfadjoint Dirac systems had been actively 
studied since 1950s [37,39] and various interesting results were published last 
years (see, e.g., [4,21-24,30,34,43,58]). In particular, the inverse problem 
to recover the m x m locally square summable potential V of the Dirac 
system (2.1) from its Weyl function (p{z) was dealt with in [54]. The case 
of rectangular (not necessarily square) matrix functions v was studied there 
and, moreover, only local square-summability of v was required, which was 
essentially less than in the preceding works. Thus a problem formulated by F. 
Gesztesy was solved and interesting applications to the case of Schrodinger- 
type operators with distributional matrix-valued potentials followed [25]. 

Theorem 2.3 [54] Let Dirac system (2.1) he given on [0, cxd), let its po¬ 
tential V be locally square-summable and let ip be the Weyl function of this 
system. Then V is uniquely recovered from ip. 

It is essential that a procedure to solve inverse problem is given in [54] (and 
will be formulated in the next subsection). Therefore, we see that if we know 
the evolution (p(t,z) of the Weyl function of the system y^ = G{x,t,z)y, 
where G is given by (2.5), we may recover the solution v{x,t) of dNLS (2.3). 
The same scheme works for many other integrable equations. 
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2.2 Recovery of the potential from the Weyl function 

In this subsection we describe a procedure to solve inverse problem for system 
(2.1) with locally square summable potentials, which is summed up in [54, 
Theorem 4.4]. The potential v of (2.1) is easily expressed in terms of the 
block rows (3 and 7 of the fundamental solution u{x, z) (normalized by (2.6)) 
at z = 0: 


/3{x) = 0]M(a;,0), 7 ( 0 ;) = [O ^(a;, 0). (2.8) 

Indeed, (2.1) yields u{x, z)*ju{x,z) = j. Hence, u{x,0)ju{x,0)* = j, and so 
I3j^* = lrm, 7J7* = -^m 2 , I3j^* = 0. (2.9) 

Using (2.1) and (2.9) we derive 

v{x) = if3'{x)j^{x)* (^f3'{x) := • (2-10) 

Moreover, relations (2.9) and initial condition /1(0) = \lmi O] (which follows 
from (2.6), (2.8)) provide a procedure to recover (3 from 7 . Namely, we have 

I3{x) = I3i{x)'i3{x), 'l3{x):=[lmi 7i(a^)*(72(a:)*)"^] ; (2.11) 

= -m' 3 n{W*r\ m) = in.,, ( 2 . 12 ) 


where Pi is the left mi x mi block of P and 71 and 72 are the m 2 x mi and 
m 2 X m 2 , respectively, blocks of 7 . We note that (2.9) yields det 72 7 ^ 0. It 
remains to recover 7 from ^p. Now, we formulate the corresponding theorem. 


Theorem 2.4 Let Dirac system (2.1) he given on [0, cxd), let its potential V 
be locally square-summable and let <p be the Weyl function of this system. 

Then V is uniquely recovered from 7 using formulas (2.2) and (2.10), 
where P is given by ( 2 . 11 ) and Pi {in ( 2 . 11 )) is the unique solution of the 
first order linear differential system {with initial condition) in (2.12). In 
turn, the matrix function 7 is uniquely recovered from <p in three steps which 
are described below. 

{i) Introduce a matrix function <hi(a;) taking a Fourier transformation 


$1 



-e’^^l.i.m.a^oo 

IT 


e-k^ 


Tii + 17 ) 

2 i(^ + i^) 




7 > 0, 


(2.13) 


—a 
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where l.i.m. stands for the entrywise limit in the norm of 0 < 

i < oo. Here $i(a;) (0 < a; < cxd) is a well-defined differentiable matrix 
function {with a locally square-summable derivative), which does not depend 
on rj > t). 

{ii) Introduce a family of pairs of operators 11; G b(^C^, and Si 

acting in (/ G M+) : 


Tlig = ^i{x)gi + g2 for gu G C™^ g = 


91 

92 


I (*X+t 


S, = I-^ 





0 J\x—t\ 


r -I X — t 




r -\-t — X 


dr ■ dt. 


(2.14) 

(2.15) 


The defined above operator Si belongs B{j^‘f^^{0,l)j, it is positive-definite 
{i.e., Si > 0 ) and boundedly invertible, and the matrix function is 

absolutely continuous with respect to 1. 

{Hi) The matrix function {Hamiltonian) H{1) = 7 (/)* 7 (/) is recovered via the 
formula 


m) = (n;sr'n,)'. 


(2.16) 


From H we recover first 


72 Si 


[O 1 ^ 2 ] H 



-1 

[0 1 ^ 2 ] H 



(2.17) 


Using 72 Si? we recover 72 from the equation 

72 = 72 ( 72 "Si)'( 72 “Si )*(^™2 - ( 7 ^Si)( 7 ^Si)*) 72 ( 0 ) = Im 2 - (2-18) 

Finally, from 7 ^^ 7 i and 72 the expression for 7 is immediate. 

We note that indices in the notations of operators (e.g., in the notations of 
operators fl; and Si above as well as in the notations S^ and Si in The¬ 
orem 3.5) indicate the spaces, in which the operators act, and should be 
distinguished from the indices meaning partial differentiation like in (2.3). 
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2.3 Zero curvature representation as a compatibility 
condition 

In this subsection we do not require any special structure of G and F (like the 
structure in (2.5)) and consider a general-type zero curvature representation 
(2.4). The term ’’compatibility condition” is very often used with respect to 
the representation (2.4) and it is easy to derive (2.4) from the existence of 
common fundamental solutions w (i.e., from the compatibility) of systems 

Wx = Gw, Wt = Fw. (2-19) 

In other words, it is easy to show that (2.4) is a necessary compatibility 
condition for systems (2.19). However, the proof of sufficiency is more com¬ 
plicated. It appeared hrst in [59, 60] and in greater detail and generality 
in [53,58]. 

More precisely, assuming that (2.4) holds in the semistrip 

= {(a^, t) : 0 < X < oo, 0 < t < a} (2.20) 

we introduce solutions of systems Wx = GW with hxed values of t and of 
systems Rt = FR with hxed values of x: 

Wx{x,t,z) = G{x,t,z)W{x,t,z), W{0,t,z) = Im] (2.21) 

Rt{x,t, z) = F{x,t, z)R{x,t, z), R{x,0,z)=lm- (2.22) 

Next, we formulate [58, Theorem 6.1]. 

Theorem 2.5 Let m x m matrix functions G and F and their derivatives 
Gt and Fx exist on the semistrip Hq, let G, Gt, and F be continuous with 
respect to x and t on fla, and let (2.4) hold. Then we have the equality 

W{x, t, z)R{0, t, z) = R{x, t, z)W{x, 0, z). (2.23) 


Thus, (2.4) implies (2.23), and so 

w{x, t, z) := W{x, t, z)R{0, t, z) = R{x, t, z)W{x, 0, z) 


(2.24) 



is the fundamental solution of systems (2.19) normalized by the initial con¬ 
dition 


w{Q,Q,z)=I^. (2.25) 

The problem of compatibility of the systems (2.19) may be formulated 
as an initial-boundary value problem. Indeed, normalizing w in (2.19) via 
(2.25) and taking into account (2.21) and (2.22) we obtain initial-boundary 
conditions on w. 

w{x, 0, z) = W{x, 0, z), tc(0, t, z) = R{0, t, z). (2.26) 

From this point of view, formula (2.24) gives a solution of the initial-boundary 
value problem (2.19), (2.26). This solution plays a crucial role in the study 
of the evolution of Weyl functions. 

2.4 Evolution of Weyl functions 

Inverse Spectral Transform (instead of the Inverse Scattering Transform) was 
hrst used for initial-boundary value problems in [15,35]. More precisely, a 
special kind of initial-boundary value problem for Toda lattice with a linear 
law of evolution of the spectral and Weyl functions was studied in [15,35]. In 
the papers [59,60], an essentially more general case of the initial-boundary 
value problem for Toda lattice as well as some initial-boundary value prob¬ 
lems for continuous integrable systems (including square matrix dNLS) were 
dealt with, and the law of evolution of the Weyl function was presented in 
the form of Mobius transformation. See further results and references, for 
instance, in [46,48,58]. 

In this section we demonstrate Inverse Spectral Transform approach on 
the important case of dNLS (2.3), where v is an rrii x m 2 matrix function. 
For that purpose, we use W and R given by (2.21) and (2.22) assuming that 
G and F have the form (2.5). First, we omit t in the notations and introduce 
the set (Weyl circle) Af{b, z) of functions of the form 

0(6, P) = [0 W{b, z)-^P{z) ( [Irn, 0] W{b, z)-^P{z)) (2.27) 
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where V{z) are m x mi nonsingular meromorphic matrix functions with 
property-j, that is, 

V{z)*V{z) > 0, V{z)*jV{z) >0 {ze C+). (2.28) 

The next formula gives an important property of Weyl functions (p{z): 

(p{z) = lim (l){b,z) (2.29) 

6—>-CXD 

for any set of functions z) G A/'(6, z). In order to derive an expression for 
the Weyl function (p(t, z) of the system = G{x, t, z)y, we rewrite (2.23) in 
the form of the equality 

W{b, t, z)-^ = R{0, t, z)W{b, 0, z)-^R{b, t, z)-^ 

and substitute this equality into (2.27). Now, passing to the limit and using 
(2.29) we obtain the required expression for (p(t,z). That is, we obtain the 
dependence of the Weyl functions on the parameter t (i.e., evolution of the 
Weyl function) in the case of dNLS. 

Theorem 2.6 [56] Let an mi x m 2 matrix function v{x,t) be continuously 
differentiable on the semistrip hla and let v^x exist. Assume that v satisfies 
the dNLS equation (2.3) as well as the following inequalities {for all 0 < t < a 
and some values M{f) G M+) : 

sup \\v{x,s)\\ < M{t). (2.30) 

0<s<t 

Then the evolution (p{t,z) of the Weyl functions of Dirac systems 
yx{x,t,z) = G{x,t, z)y{x,t, z) is given {for z G C+) by the equality 

ip{t,z) = {R 2 i{t,z) +R 22 {t,z)ip{t),z))[Rii{t,z) +Ri 2 {t,z)ip{t),z)) \ (2.31) 

where Rik{t,z) are the m* x blocks of R{0,t, z). 

3 Generalized Weyl functions 

3.1 Skew-selfadjoint Dirac system 

The system 

-^y{x,z) = {izj+jV{x))y{x,z) (x > 0, ^ G C), (3.1) 
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where j and V have the same form (2.2) as in the Dirac system (2.1), is 
called skew-selfadjoint Dirac. Weyl theory of skew-selfadjoint Dirac systems 
was studied in the papers [22,28,46,48]. Some further references as well 
as the results of this subsection are contained in [58, Ch.3]. Like in the 
case of selfadjoint Dirac system, the fundamental solution of system (3.1) is 
denoted by u{x,z), and this solution is normalized by the condition (2.6). 
The notation Cm stands for the open half-plane {z : 9^(2:) > M > 0}. 

Definition 3.1 An m 2 x mi matrix function (p, which is holomorphic in Cm 
{ for some M > 0) and satisfies the inequality 

\lmi p>{z)*^.u{x^z)*u{x^z) dx<oo, ZECm, (3.2) 

lP>[z)_ 

is called a Weyl function of the skew-selfadjoint Dirac system (3.1). 

First, assume that v is bounded, that is, 

||n(a;)|| < M for a; G [0, cxd). (3.3) 

For 2 : G Cm (with M in Cm given in (3.3)), we have the following proposition. 

Proposition 3.2 Let Dirac system (3.1) be given on [0, cxd) and assume that 
(3.3) holds. Then there is a unique Weyl function ip{z) {for z G Cm) of this 
system. Moreover, ip{z) is contractive in Cm and the inequality 

sup e~^^^u{x, z) < 00 (3.4) 

x<1,zGCm }P{^). 

holds on any finite interval [0, /]. 

Definition 3.3 A generalized Weyl function {GW-function) of the system 
(3.1), where v is locally bounded on [0, cx)), is an m 2 x mi matrix function p> 
such that for some M > 0 it is analytic in Cm and the inequality (3.4) holds 
for each I < 00 . 

Our next proposition (together with Proposition 3.2) shows that, in the case 
of V bounded on the semiaxis, the dehnitions of Weyl function and GW- 
function are equivalent, that is, GW-function coincides with Weyl function. 
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Proposition 3.4 For any system (3.1), where v is locally bounded on [0, cx)), 
there is no more than one GW-function. 

The following theorem (see [58, Theorem 3.30]) is necessary in the process 
of solving Goursat problem for the sine-Gordon equation in Section 4. The 
procedure to construct (solve inverse problem) given in this theorem 

is in many respects similar to the corresponding procedure for the selfadjoint 
Dirac system but we do not require here apriori that is a GW-function. 
We note that fd and 7 stand again for the block rows of m, that is, (2.8) holds. 

Theorem 3.5 Let an m 2 x mi matrix function ip{z) be holomorphic in Cm 
{ for some M > 0) and satisfy condition 

sup \\z^{(p{z) - 0oA)|| <00 {z e Cm), (3.5) 

where 0o is an m 2 x mi matrix. Then ip is a GW-function of a skew- 
selfadjoint Dirac system. This Dirac system is uniquely recovered from p 
using the following procedure. 

{i) First recover the matrix function <hi : 

\2J TT J_a 2 i (^ + 17 ) 

where l.i.m. stands for the entrywise limit in the norm T^( 0 , £) (0 < £ < cx). 
{ii) Next, for the values I G M+, introduce operators Si E B (L^(0,/)) : 

pi pmm(x,t) 

5'; = /+ / s{x,t) ■ dt, s{x,t) = / ^[{x — r)^[{t — r)*dr. (3.7) 

Jo Jo 

These operators are well-defined, positive definite and boundedly invertible. 
{Hi) We recover /3{x) via the formula : 

fi{x)=[lm, 0]-j^\sfi^^'i){t)*[^i{t) Im,]dt, (3.8) 

where is applied to columnwise. Since fdfd* = Imi, we can construct 
a differentiable matrix function 7 {with locally bounded derivative) such that 

= 0, 77 * > 0, 7(0) = [0 Im,] ■ (3.9) 
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Then 7 = '&'y, where h} is determined by the equation and initial eondition 
below : 


=-^7'7*(77*) \ i9{0) = 1^2- (3.10) 

{an) Finally, v is given by 

v{x) = f3\x)'y{x)*, (3.11) 

V has the form (2.2), and both v and V are locally bounded. 

3.2 Linear system auxiliary to the N-wave equation 


Nonlinear optics (iV-wave) equation 

[D, gt] - [D, = [D, g] [D, g] - [D, g] [D, p], (3.12) 

where g{x, t) = g{x, t)* is an m x m matrix function, [D, p] = Dg — gD and 

D = diagjdi, ^ 2 ,..., dm}-, di > d 2 > ■ ■ ■ > dm > 0, (3.13) 

D = D* = diag{di,d2, - - - ,dm}, (3.14) 

admits zero curvature representation (2.4), where G and F have the form 

G{x,t, z) = izD — ({x,t), F{x,t, z) = izD — ({x,t)] (3.15) 

C:=[^,^], C-=[D,g]-, (3.16) 

see [62] for the case N = 3 and [1] for N > 3. 


We shall need some preliminary results on the Weyl theory of the aux¬ 
iliary system y^ = Gy [x > 0) from [58, Ch. 4], see also [47,51]. The 
normalized fundamental solution u of such system is dehned by the formula 

-^u{x,z) = {izD - C,{x))u{x,z), u{d,z) = Im (C = -C*)- (3.17) 

Here and further we assume that is a hxed matrix satisfying (3.13). 
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Definition 3.6 A generalized Weyl funetion {GW-funetion) of system (3.17), 
where ( is locally bounded, is an mxm matrix function (p such that for some 
M > 0 it is analytic in the domain Cfj = {z : ^5(2:) < — M} and the inequal¬ 
ity 

sup z)(p{z) e'x.Y>{—izxD}^ < 00 (3.18) 

x < l , ^{ z)<-M 

holds for each I < 00 . 

The inverse spectral problem for system (3.17) is the problem to recover 
(from an analytic matrix function (p{z)) a locally bounded potential ({x) = 
—({x)* (Cfcfc = 0) such that 99 is a GlT-function of the corresponding system 
(3.17), that is, (3.18) is valid. The notation fOT stands for the solution of this 
inverse problem, that is, for an operator mapping the pair D and (p into the 
corresponding potential C (i.e., iM{D,ip) = Q. 

Theorem 3.7 For any matrix function <p{z) which is analytic and bounded 
in Cfj = {z : 7 ^( 2 ;) < —M} and has the property 

00 

J {(p{^ + 17]) - I^y {p>{^ + ip) - df < 00 {7]<-M), (3.19) 

—00 

there is at most one solution of the inverse spectral problem. 

An existence condition (for the solution of the inverse spectral problem) is 
given in the next theorem. 

Theorem 3.8 Let (p(z) be analytic in Cfj for some M > 0 and satisfy in 
the inequalities 

sup II z{ip{z) — Im) II < cxD, det p){z) 7 ^ 0. (3.20) 

Assume also that for some matrix 00 and for all fixed rj < —M we have 

{f + ip) {(p{^ + ip) -Im- 0o/(^ + ih)) e L^^^(-oo, 00 ). (3.21) 

Then the solution of the inverse problem exists and is unique. 
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Remark 3.9 The procedure to construct Tl{D, ip) under conditions of The¬ 
orem 3.8 is given in [58, Theorem f.lO], This procedure is similar to the 
procedures in Theorems 2.f and 3.5. 

Further in this paper, we deal with the case of C bounded on [0, cxd): 

sup ||C(a;)|| < oo, (3.22) 

0<x<oo 

and GhF-functions ip{z) = normalized by 

Tij{z) = l for i = j, ip.j{z)=0 for i>j. (3.23) 

The next proposition follows from [58, Subsections 4.1.1 and 4.1.3]. 

Proposition 3.10 When (3.22) holds, a normalized {by (3.23)) GW-function 
ip{z) of (3.17) exists and is unique. 


4 Sine-Gordon theory in a semistrip 


1. Sine-Gordon equation in the light cone coordinates (SGE) has the form 

^V- = 2sm(2«). (4,1) 

Local solutions of the Goursat problem for SGE were studied in [38,40] and 
global solutions were constructed hrst in [48]. The results of this section are 
obtained mostly in [48] with some modihcations and developments in [58, 
Section 6.2] (see also further references in [58]). SGE admits zero curvature 
representation (2.4), where G and F have the form [2]: 


G{x,t,z) = izj + jV{x,t), V= ^ ^ 

cos(2-^) sin(2-^) 

sin(2-^) — cos(2-^) 


F(x,t,z) = — 

iz 


^ I 

V = 

ox 




i,, (4.2) 


(4.3) 


and j has the form (2.2) with mi = m 2 = 1. Hence, system = Gy is 
a skew-selfadjoint Dirac system (3.1), where mi = m 2 = 1. The evolution 
ipit, z) of its GfF-function is given by the formula 


R 2 l{t, z) + R 22 {t, z)ip{0, z) 
Rii{t,z) -T Ri2{t,z)(p{0,z)' 


(4.4) 
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where, as usual, Rik{t, z) are the blocks of -R(0, t, z) determined by (2.22). If 
the boundary value ■^(0, t) is given, using (2.22) and (4.3) we recover Rik(t, z). 
Let us formulate the evolution result rigorously. 

Theorem 4.1 Let'ip{0,t) and be continuous functions on [0, a) and 

hla; respectively, and let tjjxtixR) exist. Let also SGE (4.1) hold on Lla, and 
assume that 99 ( 0 , z) is the GW-function of the system y^ = G{x, 0, z)y, where 
G is given by (4.2). 

Then, the function (p(t,z) of the form (4.4) is the GW-function of the 
system y^ = G{x, t, z)y. 

If the initial-boundary conditions 


'ip{x,0) = hi{x), 'ip{0,t) = h 2 {t), hi(0) =/i2(0) (4.5) 


{hk = hk for {k = 1,2)), are given, we recover (p{f),z) from the first condi¬ 
tion, recover Rik(t,z) from the second condition, and construct (p(t,z) using 
(4.4). Then we recover v{x,t) using the procedure to construct the solution 
01t(</9) of the inverse problem, which is described in Theorem 3.5. From v 
we immediately recover ip and show that ip satishes (4.1) and (4.5) (see the 
proof of [58, Theorem 6.19]). Thus, an existence theorem follows. 


Theorem 4.2 Assume that hi is boundedly differentiable on all the finite 
intervals on [0, cxd) and that ^2 is continuous on [0, a). 

Moreover, assume that the GW-function (p{f),z) of the system (3.1), 
where 


mi = m 2 = 1, V (x) 


0 h[ (x) 
h[ (x) 0 


(4.6) 


exists and satisfies (3.5). Then a solution of the initial-boundary value prob¬ 
lem (4.1), (4.5) exists and is given by the eguality 


ip{x,t) = h 2 (t) 




(4.7) 


where <p(t,z) is obtained from (4.4). 


Remark 4.3 Under conditions of Theorem f.2, for each 0 < c < a there is 
M{c) > 0 such that all functions <p(t, z) (0 < t < c) satisfy GW-function re- 
guirement (3.4) and asymptotic ineguality (3.5) in the same half-plane Cm{c)- 
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Some sufficient conditions on hi, under which the requirements on (/)(0, z) 
in Theorem 4.2 hold, were derived using important paper [11] (see [58, Corol¬ 
lary 6.21]). We formulate these conditions. 

Corollary 4.4 If v E then there is a GW-function ip of the system 

(3.1) {where V has the form (2.2) and mi = m 2 = 1). Moreover, if v is two 
times differentiable andv, v', v" G L^(M+), then this GW-funetion <p satisfies 
the asymptotic condition (3.5). 

2. Complex sine-Gordon equation was introduced (and its integrability was 
treated) [41,45] only several years after the seminal paper [2] on the integra- 
bility of SGE was published. Complex sine-Gordon equation is more general 
than SGE and has the form 

+4cosi/'(sin</l)"’xrXi = 2sin(2i/>), Xxt = (2/sin(2j/i))(i/>iX, + V’ai), 

(4.8) 

where tf = and y = y. There are also two constraint equations 

2(cos'^)^ya; — (sin'^)^a;j; = 2c(sin'^)^, 2(cos'^)^yt -|- {sinipYut = 0, (4.9) 

where c is a constant (c = c = const) and u = u. For further develop¬ 
ments and applications of the results on (4.8), (4.9) see, for instance, some 
discussions and references in [10,20,58]. Below we formulate several results 
from [58, Subsection 6.2.1]. 

If sin(2-^) 7 ^ 0 and (4.8) and (4.9) hold, then zero curvature equation 
(2.4), where 


G{x, t, z) 

:= izj + 

jV{x,t), 

F{x,t,z) 

:= ^{x,t)*j^{x,t), 

z c 

( 4 . 10 ) 

■ 1 

^=[o 

0 

-1 J ’ 

V 

V 

0 J ’ ^ ^ 


-2x) 






( 4 . 11 ) 

■d{x,t) : = 

Di{x,t) 

cosfi 

isin'^ 1 
/ F>2 

{x,t), 

( 4 . 12 ) 

1 SlU'^ 

cosip 

Di = exp 

{i(y -|- ( 

^/2))j}, 

D 2 = exp 

{i(y- (cn/2))j}. 

( 4 . 13 ) 
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also holds. Now, we introduce initial-boundary conditions: 


v{x,Q) = hi{x), =/i 2 (t), x(0,t) =/i 3 (t), a;(0,0) = /i 4 . (4.14) 


Evolution of the GhE-function 99 is given in the next theorem. 


Theorem 4.5 Let uj{x,t)} he a triple of real-valued and 

twiee continuously differentiable functions on Lla- Assume that sin( 2 -^) 7 ^ 0 , 
thatv given by (4.11) is bounded on Lla, and that complex sine-Gordon (4.8), 
(4.9) and conditions (4.14) hold. 

Then the GW-functions g:>{t,z) of the auxiliary skew-selfadjoint Dirac 
systems yx = Gy, where G{x,t,z) is defined via (4.10) and (4.11), exist in 
some Cm and have the form (4.4). 

Here R(t,z) = R{0,t,z) = [Rik(t, z)'^. is defined by the equalities 


d . 1 




cos(2/i2(t)) isin(2h2(t)) 

[ —isin(2h2(t)) — cos(2/i2(t)) J 




R{0,z) = l 2 , d{t) := hsiO) --h4 + h'fif) {sin h2{0) df] 


and 99 ( 0 , ; 2 ) is the GW-function of (3.1) where mi = m 2 = 1 and v = hi 


Corollary 4.6 There is at most one triple {fi{x,t), x(a;, t), u{x,t)} of real¬ 
valued and twice continuously differentiable functions on Da such that 
sin( 2 -^) 7 ^ 0, that v is bounded, and that complex sine-Gordon (4.8), con¬ 
straints (4.9) and initial-boundary conditions (4.14) are satisfied. 


5 Reduction of the initial-boundary condi¬ 
tions and unbounded solutions 

5.1 Reduction of the initial-boundary conditions in a 
quarter-plane 

1. In this subsection we consider integrable wave equations in the quarter- 
plane 

= {(a^, t) -. 0 < X < 00 , 0 < t < cxd}. (5.1) 
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First, we consider the defocusing NLS (2.3) with the initial-boundary condi¬ 
tions 

v{x,^) = hi{x), n(0,t) =/i2(t), = hzit). (5.2) 

Given hi{x), we recover 99 ( 0 , z) via (2.29) and given h 2 {t) and h^{t) we recover 
i?(0,t, 2 ;) (i.e., Rik(t,z) where i,k = 1,2) using (2.22). Hence, given the 
initial-boundary conditions hk {k = 1,2, 3) we recover the right-hand side of 
(2.31) and obtain the evolution ip(t, z) of the Weyl function. 

Theorem 5.1 [56] Let v satisfy the conditions of Theorem 2.6 and boundary 
conditions given by the second and third equalities in (5.2). Moreover, let the 
boundary value functions h^ and h^ be continuous and bounded, i.e., 

sup ||/i 2 (t)|| < M, sup ||/i 3 (t)|| < M (5.3) 

0<t<C3O 0<t<CxD 

for some M, M E M+. Then, in the domain 

V = {z: ^{z) > 1/2, 3fJ(;2) < -M} (5.4) 

we have the equality 

(p{fl,z) = - lim R 22 {t,z)~^R 2 i{t,z). (5.5) 

t —>-00 

Scheme of the proof. The contractiveness of ip{t, z) yields 

> 0 . (5.6) 

Putting R(t,z) := R{0,t,z) and using formula (2.31) for ip(t,z), we derive 

R{t,z) {Rn{t,z) + Ri2{t,z)(p{0,z)). (5.7) 

L<p(0,^)J 

It is immediate from (5.6) and (5.7) that 

[im, ^{0,zy]R{t,z)*jR{t,z) >0- (5-8) 
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(5.9) 


Moreover, the first inequality in (5.3) implies that 


— [ —R{t, z)*jR{t, z)) > MR{t,z)*R{t,z), z eV. 


From (5.8) and (5.9) we see that 

/ OO 

[Ira^ V^(0, Z)*] R{S, ZyR{s, z) 


.<^( 0 ,^) 


ds < (1/M)J„ 


(5.10) 


After some considerations, using the boundedness of the left-hand side of 
(5.10) and the inequalities (5.3), one can show that 


lim 

t^OO 


R{t, z) 


-‘■mi 

[(^(0,^) 


= 0 . 


(5.11) 


On the other hand, it easily follows from (5.9) that 


R 22 (t, z)*R 22 (t, z) > Im 2 - 

Therefore, formula (5.11) implies (5.5). An analogue of (5.5) for a scalar 
dNLS (and with a somewhat different proof) appeared already in [59]. 

In view of analyticity of Weyl functions, formula (5.5) means that 99 ( 0 , z), 
ip(t, z), and thus also v(x, t), may be recovered from the boundary conditions. 
2. A similar result is valid for the focusing nonlinear Schrodinger equation: 


2vt + i(vxx + 2vv*v) = 0, (5-12) 

where v is an mi x m 2 matrix function. Equation (5.12) admits [64, 65] 
representation (2.4) where 

G = izj+jV, F = i{z^j-izjV-{Vx+jV^)/2), (5.13) 

and j and V are dehned in (2.2). Evolution of the Weyl function is described 
in this case by the following theorem. 

Theorem 5.2 [28] Let an mi x m 2 matrix function v{x,t) be continuously 
differentiable on hla and let v^x exist. Assume that v satisfies the focusing 
nonlinear Schrodinger eguation (5.12) as well as the following inegualities : 

sup ||t;(a;,t)|| < M, sup ||v 2 ;(a;,t)|| < cx) foreach 0 < c < a. 
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Then, the evolution (p(t, z) of the GW-funetions of the skew-selfadjoint Dirae 
systems = Gy is given {for z G Cm) via the Mobius transformation (2.31), 
where the eoeffieients Rik are determined by the formula (2.22) and the equal¬ 
ity {Rik{t, z)}"^:= R{0,t,z). Here G and F have the form (5.13). 

Remark 5.3 It is immediate that if the conditions of Theorem 5.2 hold for 
each a < oo, then (2.4) describes the evolution of (p{t, z) for allO <t < oo. 

In order to derive an analogue of Theorem 5.1 we shall require boundedness 
of V on hloo (instead of boundedness on all fla) and boundedness of ^2,(0,^) 
on M+: 

sup ||n(a;, f)|| < M, sup ||n2;(0,t)|| < cxd. (5-14) 

(a;,t)GOoo iGK+ 

Then, using (2.22) (see, e.g., [28, formula (4.14)]), we derive inequality (5.9) 
for R{t, z) = R(0, t, z) in a domain 

V = {z: Q(z) > Ml, 3fJ(;2) > M 2 }, (5.15) 

and with some M > 0, Mi > M, M 2 > 0. After that, the same considerations 
as in the scheme of the proof of Theorem 5.1 yield our next theorem. 

Theorem 5.4 Let v satisfy the conditions of Theorem 5.2 {for each a G M+) 
and inequalities (5.14). Then, in some domain (5.15) we have the equality 
(p{0,z) = -\imt-,ocR 22 {t,z)-^R 2 i{t,z). 

3. Finally, let us consider SGE (4.1), where v = is bounded in Goo: 

sup \'il)x{x,t)\ < M. (5.16) 

(x,t)efioo 

The next proposition follows (see [58, Corollary 6.25]) from Theorem 4.1. 

Proposition 5.5 Let fj given in Coo satisfy inequality (5.16). Assume also 
that the conditions of Theorem 4-1 and initial-boundary conditions (4.5) hold 
for each (a G M+). Then, for values of z E Cm, such that the inequalities 

(cos(2h2(t)) — e(2:))A(2:) > |3fJ(2:) sin(2/i2(t))| (5-17) 

are valid for some e{z) > 0 and for all t >0, we have 

99 ( 0 , z) = - lim R 2 i{t, z)/R 22 {t, z). (5.18) 

t—>-oo 
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5.2 Unbounded solutions of SGE 

Theorem 4.2, Corollary 4.4 and Proposition 5.5 allow us to construct wide 
classes of unbounded (in hloo) solutions of SGE, that is, solutions which 
do not satisfy (5.16). The simplest case is the case h- 2 {t) = 0 (see [58, 
Proposition 6.27]). In this case, (5.17) holds for all 2 : G Cm {M > 0) and, 
if other conditions of Proposition 5.5 are valid, (5.18) yields v?(0,2:) = 0. 
Since 9H(0) = 0, we derive hi = 0. In other words, any solution of SGE 
with /i 2 (t) = 0 and hi{x) ^ 0 does not satisfy conditions of Proposition 5.5. 
Solutions of initial-boundary value problems for SGE may be constructed 
using Theorem 4.2 and Corollary 4.4. 

Proposition 5.6 Assume that hi{x) = hi{x) ^ 0 is three times differen¬ 
tiable for X >t), that 

h[, hf hf e L\R+), hi(0) = 0, 

and that /i 2 = 0. Then one can use the procedure given in Theorem f.2 in or¬ 
der to construct a solution f of the initial-boundary value problem (4.1), (4.5) 
for SGE, and the absolute value of the derivative {for the constructed ip) 
is always unbounded in the quarter-plane hloo- 

Indeed, the solution constructed in Proposition 5.6 satishes the conditions of 
Proposition 5.5 excluding, possibly, (5.16), and so (5.16) does not hold. 

Some classes of unbounded solutions of the KdV equation with the mi¬ 
nus sign before the dispersion term are constructed in [52] using low energy 
asymptotics of the Weyl functions. 

5.3 Reduction of the initial-boundary conditions in a 
semistrip 

1. In this paragraph we present some results from [55, Section 4] on the 
TV-wave equation (3.12) in the semistrip hla- By ip{t,z) we denote evolution 
of the GhP-function corresponding to this equation or, more precisely, (p{t, z) 
are GhP-functions of the systems yx = G{x,t,z)y, where G is expressed 
via g using (3.15) and (3.16). The matrix function R{t,z) = R{0,t,z) is 
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determined via (2.22), where F (similar to G) is given by (3.15) and (3.16). 
First, we express normalized (by (3.23)) GhF-functions (p(t,z) via 9?(0,;2) 
and R{t, z) (in other words, we derive the evolution of the normalized GW- 
functions). Next, we consider the subcase when the entries of D are ordered 
in the same way as the entries of D: 

D = diagjdi, ^2, • • •, dm}-, di > ^2 > • • • > > 0. (5.19) 

Theorem 5.7 Let g = g* satisfy the N-wave equation (3.12), where D and 
D have the form (3.13) and (3.14), respectively. Assume that g{x,t) is uni¬ 
formly bounded and continuously differentiable on 
Then the matrix functions 

ijk{t,z) :=[lk t)\R{t,z)y:>{f),z) 

■^m—k 

x([0 I^_k]R{t,zM0,z)\ ]] , (5.20) 

\ L -'m-fc J / 

where v^(0,2:) is the normalized GW-function of the system y^ = G{x,0,z), 
are well-defined for 1 < k < m. The normalized GW-functions (p(t,z) are 
given {in Cfj for some M > 0) by the formula 

’ 1 ' 

{^i^k+i{t,z)}^^^^=fik{t,z) {zeCf^), (5.21) 

0 

and by the normalization conditions (3.23). 

From Theorems 3.7 and 5.7 follows the next result. 

Theorem 5.8 For the case where the entries of the matrix D in (3.12) 
are positive and ordered as in (5.19), there is no more than one uniformly 
bounded and continuously differentiable on hla solution g = g* {of the N- 
wave equation (3.12)), having the initial values p(a;, 0) such that (p{t),z) is 
bounded and (3.19) holds. That is, there is no more than one solution of the 
corresponding initial value problem. 
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Taking into account the results on fundamental solutions from [11], the re¬ 
quirements on 99(0, z) given in Theorem 5.8 are reformulated below in terms 
of the sufficient requirements on the initial condition 

f?(a;,0) = p{x) = p{x)* {0 < X < 00 ). (5.22) 

Proposition 5.9 Suppose that the initial condition p{x) is absolutely con¬ 
tinuous on [0, cxd) and the entries of p{x) and p'{x) belong Then, 

the normalized GW-function 99(0,2;) of the system 

y'{x,z) = {izD - C,{x))y{x,z) (x > 0), C=[D,p] (5.23) 

is analytic and bounded in {for some M > 0) and (3.19) is valid. 

Theorem 5.8 is proved in [55, Section 4] using the fact that 

sup \\R{t, z)p{0, z) exp{—iztD}\\ < 00 , (5-24) 

tG[0,a), ^(z)<—M 

which yields (see [58, pp. 108, 109]) that 99(0,2:) determines not only the 
initial but also the boundary values of g. In particular, when 99(0, z) satishes 
the conditions of Theorem 3.8, then, according to Remark 3.9, we have a 
procedure to recover the initial and boundary conditions from 99(0, z). More 
precisely, we have 

ax,0) = m{D,p{0,z)), ao,t) = m{D,<p{o,z)), (5.25) 

and, using (3.16), we easily obtain ^?(a;, 0) and ^?(0,^). Indeed, without loss 
of generality we may assume that all the entries of g on the main diagonal 
are equal to zero and recover g = {gik}Tk=i C or ^ via formulas 

gu = 0, gik = Cik/{di dk') = Cik/{di dfc) for i ^ k. 

2. In this and the next paragraph we deal with the cases of dNLS (2.3) 
with quasi-analytic boundary or initial, respectively, conditions, (i.e., the 
corresponding boundary or initial value functions belong to quasi-analytic 
classes). The dehnition of quasi-analytic classes C([0,£); M) is given in Ap¬ 
pendix A. 
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Further we present some results from [55, Section 3]. Let us consider mi x 
m 2 matrix functions v{x,t), which are continuously differentiable and such 
that Vxx exists on the semi-strip hla. Moreover, we require some smoothness 
of V in the neighborhood of the point (0, 0). More precisely, we require that 
for each k there is a value £k = £k{v) > 0 such that v is k times continuously 
differentiable with respect to x in the square 

= {(a:, t) : 0 <x < Ek, 0 <t < Ek}, ^{Ek) C fia- (5.26) 

The class of such functions v{x,t) is denoted by Ce{k}a)- 

Proposition 5.10 Assume thatv{x,t) G C'£(f^a) satisfies the dNLS equation 
(2.3) on kla and that a matrix function n(0,t) orvxififi) is quasi-analytic. 

Then a matrix function v{0,t) orvx{0,t), respectively, is uniquely deter¬ 
mined by the initial condition 

v{x,t)) = h{x). (5.27) 

Proposition 5.10 is proved by presenting formulas to recover the derivatives 
(^n) (0, 0) and (^■§^Vx^ (0, 0) {k > 0) from h{x) (see the proof of [55, 
Proposition 3.2]). 

Using Proposition 2.2, Theorems 2.4 and 2.6 and Proposition 5.10 we 
derive the next theorem. 

Theorem 5.11 Assume that v G Ce(12a) satisfies the dNLS equation (2.3) 
onkla, that (2.30) holds and that the functions vif), f) and Vx{0,t) {boundary 
values) belong to some quasi-analytic classes C{[0, a); M) andC([0,a); M+), 
respectively. Then v is uniquely determined by the initial condition (5.27). 

3. Using Proposition 2.2 and Theorems 2.4 and 2.6 one can also derive 
Theorem 4.3 from [56], which we formulate in this paragraph. The proof is 
based on the formulas to recover the derivatives (0, 0) {k > 0) from 

the boundary conditions 

n(0, t) = /io(t), Vx{^,f) = hi{t) (0 < t < a). (5.28) 

Theorem 5.12 Assume that v G C£(r2a) satisfies the dNLS equation (2.3) 
on 12a, that (2.30) holds and that the initial value function v{x, 0) belongs to 
some quasi-analytic class C([0,cxd); M). Then v is uniquely determined by 
the boundary conditions (5.28). 
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6 


Dynamical Dirac system and response func¬ 
tion 


6.1 Introduction 

In this section we give several statements from our paper [57], which appeared 
in arXiv in 2015. Classical Dirac systems (2.1) are also called spectral Dirac 
systems and various new results on inverse problems for these systems ap¬ 
peared quite recently (see some references in Subsection 2.1). At the same 
time, a great and growing interest in dynamical systems and control meth¬ 
ods is reflected in the active study of the dynamical inverse problems and, 
in particular, in the study of the inverse problems for dynamical Schrodinger 
and Dirac systems [7,8,13,14,33,42] (see also the references therein). 

Dynamical Dirac system (Dirac system in the time-domain setup) was 
studied in the important recent paper [14]. The dynamical Dirac system 
considered in [14] is an evolution system of hyperbolic type and has the 
following form: 


iYt JYx + VY = 0 (x > 0, t > 0); 


Y = 

Yf 

, J = 

' 0 1' 

, V = 

'p q ' 


Y2_ 


-1 0_ 


_q -p_ 


dY 


( 6 . 1 ) 

( 6 . 2 ) 


where p = p{x) and q = q{x) are real-valued functions of x, and initial¬ 
boundary conditions are given by the equalities 


Y{x,0) = 0, a; > 0; Yi(0,f) = /(t), t>0. (6.3) 


Here / is a complex-valued function (called boundary control by [14]) and 
the input-output map (response operator) R : Yi(0, •) —)■ l2(0, •) is of the 
convolution form Rf = if r * f. Thus, R and r denote in this section 
the response operator and response function, respectively. {Y stands in this 
section for the Fourier transformation of Y, see (6.16).) 

The inverse problem consists in recovery of the potential V from the re¬ 
sponse function r. This inverse problem was considered in [14] using boundary 
control methods. 
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We note that recent results on dynamical Schrodinger and Dirac equa¬ 
tions are based on several earlier works. In his paper [17] from 1971, A.S. Bla- 
govescenskii considered dynamical system 

Yu - W. + Q{x)Y^ = 0 (6.4) 

with boundary control D(0,t) = /(t), and solved inverse problem to recover 
Q from /. A.S. Blagovescenskii established important connections between 
his problem and spectral theory of string equations. This work was developed 
further in [12] (see also references therein), where response operator appears 
in inverse problem. Finally, the inverse problem to recover the matrix po¬ 
tential Q{x) of the dynamical Schrodinger equation 

Yu - W. + Q{x)Y = 0 (6.5) 

from the response function was considered in [6] by S. Avdonin, M. Belishev, 
and S. Ivanov. 

Similar to the case of the spectral Dirac and Schrodinger equations, 
the dynamical Dirac equation is a more general object than the dynamical 
Schrodinger equation. More precisely, setting in (6.1) 

p{x) = 0, q{x) = gx{x)/g{x), where gxx{x) = Q{x)g{x), (6.6) 

and rewriting (6.1), (6.2) in the form 

(W)* = i{{Y,)x + gl2), {Y,), = i( - (W). -f glh), (6.7) 

we obtain a dynamical Schrodinger equation 

(hj)tt = {Yi)xx - {qx + q^)Y, = (W),, - {gxx/g)Yi 

= {Y,)xx-QY,. ( 6 . 8 ) 

For interesting applications of the interconnections between spectral Dirac 
and Schrodinger equations see, for instance, the papers [18,23-25,30,31] and 
references therein. 


27 



6.2 Preliminaries and estimates 

According to [14, Theorem 1], in the case where p, q, f (in (6.2) and in 
(6.3)) are continuously differentiable and /(O) = /'(O) = 0, there is a unique 
classical solution Y of (6.1), (6.3) and this solution admits representation 


Y = Yf + wf- 

y/{x,t) = fit- x) j , /(t) = 0 (t<0); 

(6.9) 

W^{x, t) = 

J X 

f{t — s)k{x, s)ds {t > X >0), 

(6.10) 

w^{x, t) = 0 

{x > t > 0), 

(6.11) 


where k{x, s) {x < s) is continuously differentiable. In particular, formulas 
(6.9) and (6.11) yield: 

Y{x,t) = 0 for 0 <t < X (hniteness of the domain of influence). (6.12) 

Representation (6.9)-(6.11) is proved in [14] using Duhamel formula. 

Let us also assume that V, / and /' are bounded: 

sup ||V(a;)|| < Ml, sup f{t) ^ < Cq, sup f{t) ^ < Cq. (6.13) 

a ;>0 t >0 _ 1 . ^>0 -L 

The following estimates are proved in [57, Section 2]. 

Proposition 6.1 Let p, q, f be continuously differentiable and let equalities 

fit) = fit) = 0 hold for t < 0. Assume that (6.13) is valid. Then the 

solution Y of the dynamical Dirac system (6.1), such that (6.3) and (6.12) 
are valid, satisfies the following inequalities 

||F(a;,t)|| < coe^*, \\Yt{x,t)\\ <CQe^\ ||Tr(a:, t)|| < Mae^*, (6.14) 

where a; > 0, t > 0, M 2 > Q is some constant, and M = 2\/2Mi. 

The kernel k of the integral operator in (6.10) satifies the inequality 


K{x,f)\\ < Me^h 


(6.15) 



In view of (6.14), we can apply to Y the transformation: 

poo 

Y{x,z)= / e^^^Y{x,t)dt, z^Cm, (6.16) 

Jo 

and Y stands in this section for the Fourier transformation of Y (which is 
taken, for the sake of convenience, for the hxed values x,z). Moreover, the 
same transformation can be applied to Yt and F^, and we have 

poo 

i / e^^^Yt{x,t)dt = zY{x, z), (6.17) 

Jo 

e^^^YJx,t)dt =—Y(x, z) =: Y\x, z), zeCm- (6.18) 

dx 

Now, applying the Fourier transformation to the dynamical Dirac system 
(6.1), we derive 

zY{x, z) + JY'{x, z) + V{x)Y (x, z) = 0. (6.19) 

Note that, according to [14], the response funetion r is given by 

r(t) = K 2 ( 0 ,t). (6.20) 

Indeed, for r of the form (6.20), using (6.3), (6.9) and (6.10) we obtain 
Fi(0,t) = fit) and 

^ 2 ( 0 , t) = i/(t) + f r{t — s)f{s)ds = iYi{0,t)+ j r(t — s)Fi(0, s)(is. (6.21) 
Jo Jo 

6.3 Response and Weyl functions 

In this subsection we always assume that the conditions of Proposition 6.1 
are fulhlled. However, it seems possible (and would be interesting) to mod¬ 
ify representation (6.9)-(6.11) for the case of matrix functions p and q and 
for weaker smoothness conditions, in which case the solution of the inverse 
problem in a much more general situation will follow. 

The equivalence transformation between spectral Dirac systems in the 
form (6.19), where J and V are given in (6.2), and in the form (2.1), (2.2) 
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(mi = m 2 = 1) is given by the relations 


y = KY, V = iq — p, K 


1 r i 1 
V 2 [-i 1, 


( 6 . 22 ) 


Since, according to [57, Section 3], u{x,z) {z G Cm) is a Weyl solution 
of (6.19) (i.e., u G /.^(O, cxd)), we see that y = KY is a Weyl solution of 

(2.1) , (2.2), where mi = m 2 = 1. Hence, in view of (6.21), dehnition (2.7) 
yields [57, Proposition 3.1]: 

Proposition 6.2 The response function r(t) of the dynamical Dirac system 

(6.1) is connected with the Weyl function (p(z) of the corresponding spectral 
Dirac system (2.1) [where mi = m 2 = 1) via the equality 


(p[z) = r[z)/[r[z) + 2i), ^ G Cm- (6.23) 


Remark 6.3 Proposition 6.2 jointly with Theorem 2.4 and the equality v = 
iq — p in the equivalence transformation (6.22) [which yields p = —iR.[v), 
q = 7y(n)) provide a procedure to solve the inverse problem of the recovery 
of the potential V of a dynamical Dirac system from its response function. 


Since mi = m 2 , it is possible to introduce a slightly different class of Weyl 
functions (that is, Weyl functions pn) via the inequality: 



i(pH[z)* Qu[x, z)*u[x, z)Q* 


0 : = 


1 

72 


h 

Ik 



z G C+, 


h 

-Wh[z) 


dx < 00 , 


(6.24) 

(6.25) 


where k = mi = m 2 . Comparing (6.24) with (2.7) we see that pu is a 
linear fractional transformation of <p. It is convenient that Weyl functions 
ifH belong to Herglotz class (i.e., i[(pH[zY — ph[z)) > 0). In our case we 
have k = 1 and the connection between (pn and r is simpler than (6.23). 
Namely, we have 


Th[z) = r[z) + i. 


(6.26) 
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Remark 6.4 Clearly, we can recovery from (pn after an easy transformation 
which maps pn into p. However, there is also an independent procedure 
in [ 49 ] {see [57, 58] as welt) to recover v from pu- Instead of the structured 
operators Si given by (2.15) and (2.13), convolution operators Si : 

~ d f’' 

S'; = — / s{x — t) ■ dt, s{x) = —s{—x)* {x > 0), (6.27) 

dx Jq 

^ J + iT])-^pH(f + ir])df^ (6.28) 

are used for this purpose. Formulas (6.26) and (6.28) imply that r{t) = 
2is'(t). Moreover, it is shown in [57] that the response function r coincides 
with a so called A-amplitude from [32] {more precisely, with the analogue of 
A-amplitude for the corresponding spectral Dirac system). We note that in 
M. G. Krein’s terminology 's' is called the accelerant. 

Interconnections between response functions, Weyl functions and A-amplitudes 
for Schrodinger equations are discussed in the interesting paper [8]. 

Some explicit solutions of the inverse problem for dynamical Dirac system 
are obtained in [57, Section 5]: 

Theorem 6.5 ]57] Let r{t) he the response function of a dynamical Dirac 
system and assume that r{f) admits representation 

r{t) = 

where the n x n (n G N) matrix a and the column vectors dj G (f = 1, 2) 
satisfy the identity 

a - a* = -i(di + ^ 2 ) + 112 )*■ 

Then the potential V of this dynamical Dirac system is given {in terms of a 
and dj) by the egualities 


V = 


p q 
9 -P. 


, p = -3f^(^'), q = ^{v)] 

v{x) = S{x)-h'^^92, A:=a + i9i{A + ^ 2 y] 

pX 

S{x) = In+ / A{t)A{t)*dt, A(t) = [e“‘*^'di e‘*^'d 2 ] • 
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A Quasi-analytic functions and matrix func¬ 
tions 

The class consists of all inhnitely differentiable on [Ojf') scalar 

fnnctions / snch that for some c(/) > 0 and for hxed constants > 0 
(/c > 0) we have 

< (c(/))^+^Mfc forall te[0,i), 0 < i < oo. (A.l) 

Recall that C(^{Mk}) is called qnasi-analytic if for the fnnctions / from this 
class and for any 0 < t < £ the eqnalities ^(t) = 0 {k > 0) yield / = 0. 
According to the famous Denjoy-Carleman theorem, the equality 

^ — = oo, Ln:= inf M\J'" (A.2) 

L/y. k>n 

n=l 

implies that the class C[{Mk}) is qnasi-analytic. 

In the case of matrix fnnctions (j)(t) (e.g, in the case of 0(t) = u(0,t)) we 
say that 0 is qnasi-analytic if the entries of 0 are qnasi-analytic and we say 
that 0 G C{\^,£)] M), where M = {Mk{i,j)}, if the entries 0*^ of 0 belong 
to qnasi-analytic classes C{{Mk{i, j)}). 

The recovery of a fnnction belonging to some qnasi-analytic class C 
from its Taylor coefficients (and from the seqnence {Mfc}) is discussed in im¬ 
portant works [9,36] and [16, Section III.8] bnt many interesting problems 
are still open. 
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